Accelerating Multiple Sequence Alignment with the Cell BE Processor

The Cell Broadband Engine (BE) Architecture is a new heterogeneous multi-core architecture targeted at compute-intensive workloads. The architecture of the Cell BE has several features that are unique in high-performance general-purpose processors, most notably the extensive support for vectorization, scratch pad memories and explicit programming of direct memory accesses (DMAs) and mailbox communication. While these features strongly increase programming complexity, it is generally claimed that significant speedups can be obtained by using Cell BE processors. This paper presents our experiences with using the Cell BE architecture to accelerate Clustal W, a bio-informatics program for multiple sequence alignment. We report on how we apply the unique features of the Cell BE to Clustal W and how important each is in obtaining high performance. By making extensive use of vectorization and by parallelizing the application across all cores, we demonstrate a speedup of 24.4 times when using 16 synergistic processor units on a QS21 Cell Blade compared to single-thread execution on the power processing unit. As the Cell BE exploits a large number of slim cores, our highly optimized implementation is just 3.8 times faster than a 3-thread version running on an Intel Core2 Duo, as the latter processor exploits a small number of fat cores.

New insights into sequence variation in the IGS region of 21 cyathostomin species and the implication for molecular identification

Cyathostomins comprise a group of 50 species of parasitic nematodes that infect equids. Ribosomal DNA sequences, in particular the intergenic spacer (IGS) region, have been utilized via several methodologies to identify pre-parasitic stages of the commonest species that affect horses. These methods rely on the availability of accurate sequence information for each species, as well as detailed knowledge of the levels of intra- and inter-specific variation. Here, the IGS DNA region was amplified and sequenced from 10 cyathostomin species for which sequence was not previously available. Also, additional IGS DNA sequences were generated from individual worms of 8 species already studied. Comparative analysis of these sequences revealed a greater range of intra-specific variation than previously reported (up to 23%); whilst the level of inter-specific variation (3-62%) was similar to that identified in earlier studies. The reverse line blot (RLB) method has been used to exploit the cyathostomin IGS DNA region for species identification. Here, we report validation of novel and existing DNA probes for identification of cyathostomins using this method and highlight their application in differentiating life-cycle stages such as third-stage larvae that cannot be identified to species by morphological means.

Bioinformatics inspired algorithm for stereo correspondence

In this paper, we exploit the analogy between protein sequence alignment and image pair correspondence to design a bioinformatics-inspired framework for stereo matching based on dynamic programming. This approach also led to the creation of a meaningfulness graph, which helps to predict matching validity according to image overlap and pixel similarity. Finally, we propose an automatic procedure to estimate automatically all matching parameters. This work is evaluated qualitatively and quantitatively using a standard benchmarking dataset and by conducting stereo matching experiments between images captured at different resolutions. Results confirm the validity of the computer vision/bioinformatics analogy to develop a versatile and accurate low complexity stereo matching algorithm.

Purification, characterization and molecular cloning of chymotrypsin inhibitor peptides from the venom of Burmese Daboia russelii siamensis

One novel Kunitz BPTI-like peptide designated as BBPTI-1, with chymotrypsin inhibitory activity was identified from the venom of Burmese Daboia russelii siamensis. It was purified by three steps of chromatography including gel filtration, cation exchange and reversed phase. A partial N-terminal sequence of BBPTI-1, HDRPKFCYLPADPGECLAHMRSF was obtained by automated Edman degradation and a Ki value of 4.77. nM determined. Cloning of BBPTI-1 including the open reading frame and 3' untranslated region was achieved from cDNA libraries derived from lyophilized venom using a 3' RACE strategy. In addition a cDNA sequence, designated as BBPTI-5, was also obtained. Alignment of cDNA sequences showed that BBPTI-5 exhibited an identical sequence to BBPTI-1 cDNA except for an eight nucleotide deletion in the open reading frame. Gene variations that represented deletions in the BBPTI-5 cDNA resulted in a novel protease inhibitor analog. Amino acid sequence alignment revealed that deduced peptides derived from cloning of their respective precursor cDNAs from libraries showed high similarity and homology with other Kunitz BPTI proteinase inhibitors. BBPTI-1 and BBPTI-5 consist of 60 and 66 amino acid residues respectively, including six conserved cysteine residues. As these peptides have been reported to have influence on the processes of coagulation, fibrinolysis and inflammation, their potential application in biomedical contexts warrants further investigation. © 2013 Elsevier Inc.

Trait directed de novo population transcriptome dissects genetic regulation of a balanced polymorphism in phosphorus nutrition/arsenate tolerance in a wild grass Holcus lanatus

The aim of this study was to characterize the transcriptome of a balanced polymorphism, under the regulation of a single gene, for phosphate fertilizer responsiveness/arsenate toler- ance in wild grass Holcus lanatus genotypes screened from the same habitat.

De novo transcriptome sequencing, RNAseq (RNA sequencing) and single nucleotide poly- morphism (SNP) calling were conducted on RNA extracted from H.lanatus. Roche 454 sequencing data were assembled into c. 22 000 isotigs, and paired-end Illumina reads for phosphorus-starved (P) and phosphorus-treated (P+) genovars of tolerant (T) and nontoler- ant (N) phenotypes were mapped to this reference transcriptome.

Heatmaps of the gene expression data showed strong clustering of each P+/P treated genovar, as well as clustering by N/T phenotype. Statistical analysis identified 87 isotigs to be significantly differentially expressed between N and T phenotypes and 258 between P+ and P treated plants. SNPs and transcript expression that systematically differed between N and T phenotypes had regulatory function, namely proteases, kinases and ribonuclear RNA- binding protein and transposable elements.

A single gene for arsenate tolerance led to distinct phenotype transcriptomes and SNP pro- files, with large differences in upstream post-translational and post-transcriptional regulatory genes rather than in genes directly involved in P nutrition transport and metabolism per se.

Cytosolic phospholipase A2 gene expression in rat mesangial cells is regulated post-transcriptionally

Cytosolic phospholipase A2 (cPLA2) is thought to be the rate-limiting enzyme in the arachidonic acid/eicosanoid cascade. The ability of various agonists to increase steady-state cPLA2 mRNA levels has previously been reported. The current study delineates the contributions of transcriptional and post-transcriptional processes to the regulation of cPLA2 gene expression in response to a variety of agonists in cultured rat glomerular mesangial cells. Epidermal growth factor, platelet-derived growth factor, serum and phorbol myristate acetate all increase the half-life of cPLA2 mRNA transcripts, indicating a role for post-transcriptional modulation of gene expression. The presence of three ATTTA motifs in the 3' untranslated region (3'UTR) of the rat cPLA2 cDNA is ascertained. Heterologous expression of chimeric constructs with different 3'UTRs ligated into the 3' end of the luciferase coding region reveals that the presence of the cPLA2 3'UTR results in reduced luciferase activity compared with constructs without the cPLA2 3'UTR. Furthermore, the luciferase activity in the constructs with the cPLA2 3'UTR is increased in response to the same agonists which stabilize endogenous cPLA2 mRNA. A negligible effect of these agonists on transcriptional control of cPLA2 is evident using promoter-reporter constructs expressed in transient and stable transfectants. Taken together, these results indicate predominant post-transcriptional regulation of cPLA2 mRNA levels.

Isolation of promoter for cytosolic phospholipase A2 (cPLA2)

Cytosolic phospholipase A2 (cPLA2) releases arachidonic acid from membrane phospholipids and is believed to be the rate-limiting enzyme in the arachidonic acid pathway. We report herein the isolation of a 3 kb fragment of rodent genomic DNA containing part of the first intron, the first exon and 5'-flanking sequence. The start site of transcription was mapped by 5'-rapid amplification of cDNA ends and corroborated by ribonuclease protection assay. The gene has a TATAless promoter with no classical Sp1 binding sites or initiator element. A microsatellite series of CA repeats was noted in the 5'-flanking region of both the rodent and human promoters. Deletion constructs have been analysed for luciferase activity and confirmed promoter activity.

Norspermidine Is Not a Self-Produced Trigger for Biofilm Disassembly

Formation of Bacillus subtilis biofilms, consisting of cells encapsulated within an extracellular matrix of exopolysaccharide and protein, requires the polyamine spermidine. A recent study reported that (1) related polyamine norspermidine is synthesized by B. subtilis using the equivalent of the Vibrio cholerae biosynthetic pathway, (2) exogenous norspermidine at 25 μM prevents B. subtilis biofilm formation, (3) endogenous norspermidine is present in biofilms at 50-80 μM, and (4) norspermidine prevents biofilm formation by condensing biofilm exopolysaccharide. In contrast, we find that, at concentrations up to 200 μM, exogenous norspermidine promotes biofilm formation. We find that norspermidine is absent in wild-type B. subtilis biofilms at all stages, and higher concentrations of exogenous norspermidine eventually inhibit planktonic growth and biofilm formation in an exopolysaccharide-independent manner. Moreover, orthologs of the V. cholerae norspermidine biosynthetic pathway are absent from B. subtilis, confirming that norspermidine is not physiologically relevant to biofilm function in this species.

The structure of an unconventional HD-GYP protein from Bdellovibrio reveals the roles of conserved residues in this class of cyclic-di-GMP phosphodiesterases

UNLABELLED: Cyclic-di-GMP is a near-ubiquitous bacterial second messenger that is important in localized signal transmission during the control of various processes, including virulence and switching between planktonic and biofilm-based lifestyles. Cyclic-di-GMP is synthesized by GGDEF diguanylate cyclases and hydrolyzed by EAL or HD-GYP phosphodiesterases, with each functional domain often appended to distinct sensory modules. HD-GYP domain proteins have resisted structural analysis, but here we present the first structural representative of this family (1.28 Å), obtained using the unusual Bd1817 HD-GYP protein from the predatory bacterium Bdellovibrio bacteriovorus. Bd1817 lacks the active-site tyrosine present in most HD-GYP family members yet remains an excellent model of their features, sharing 48% sequence similarity with the archetype RpfG. The protein structure is highly modular and thus provides a basis for delineating domain boundaries in other stimulus-dependent homologues. Conserved residues in the HD-GYP family cluster around a binuclear metal center, which is observed complexed to a molecule of phosphate, providing information on the mode of hydroxide ion attack on substrate. The fold and active site of the HD-GYP domain are different from those of EAL proteins, and restricted access to the active-site cleft is indicative of a different mode of activity regulation. The region encompassing the GYP motif has a novel conformation and is surface exposed and available for complexation with binding partners, including GGDEF proteins.

IMPORTANCE: It is becoming apparent that many bacteria use the signaling molecule cyclic-di-GMP to regulate a variety of processes, most notably, transitions between motility and sessility. Importantly, this regulation is central to several traits implicated in chronic disease (adhesion, biofilm formation, and virulence gene expression). The mechanisms of cyclic-di-GMP synthesis via GGDEF enzymes and hydrolysis via EAL enzymes have been suggested by the analysis of several crystal structures, but no information has been available to date for the unrelated HD-GYP class of hydrolases. Here we present the multidomain structure of an unusual member of the HD-GYP family from the predatory bacterium Bdellovibrio bacteriovorus and detail the features that distinguish it from the wider structural family of general HD fold hydrolases. The structure reveals how a binuclear iron center is formed from several conserved residues and provides a basis for understanding HD-GYP family sequence requirements for c-di-GMP hydrolysis.

Three motAB stator gene products in Bdellovibrio bacteriovorus contribute to motility of a single flagellum during predatory and prey-independent growth

The predatory bacterium Bdellovibrio bacteriovorus uses flagellar motility to locate regions rich in Gram-negative prey bacteria, colliding and attaching to prey and then ceasing flagellar motility. Prey are then invaded to form a "bdelloplast" in a type IV pilus-dependent process, and prey contents are digested, allowing Bdellovibrio growth and septation. After septation, Bdellovibrio flagellar motility resumes inside the prey bdelloplast prior to its lysis and escape of Bdellovibrio progeny. Bdellovibrio can also grow slowly outside prey as long flagellate host-independent (HI) cells, cultured on peptone-rich media. The B. bacteriovorus HD100 genome encodes three pairs of MotAB flagellar motor proteins, each of which could potentially form an inner membrane ion channel, interact with the FliG flagellar rotor ring, and produce flagellar rotation. In 2004, Flannagan and coworkers (R. S. Flannagan, M. A. Valvano, and S. F. Koval, Microbiology 150:649-656, 2004) used antisense RNA and green fluorescent protein (GFP) expression to downregulate a single Bdellovibrio motA gene and reported slowed release from the bdelloplast and altered motility of the progeny. Here we inactivated each pair of motAB genes and found that each pair contributes to motility, both predatorily, inside the bdelloplast and during HI growth; however, each pair was dispensable, and deletion of no pair abolished motility totally. Driving-ion studies with phenamil, carbonyl cyanide m-chlorophenylhydrazone (CCCP), and different pH and sodium conditions indicated that all Mot pairs are proton driven, although the sequence similarities of each Mot pair suggests that some may originate from halophilic species. Thus, Bdellovibrio is a "dedicated motorist," retaining and expressing three pairs of mot genes.

Roles of multiple flagellins in flagellar formation and flagellar growth post bdelloplast lysis in Bdellovibrio bacteriovorus

Bdellovibrio bacteriovorus cells have a single polar flagellum whose helical pitch and diameter characteristically change near the midpoint, resulting in a tapered wave. There are six flagellin genes in the genome: fliC1 to fliC6. Accordingly, the flagellar filament is composed of several similar flagellin species. We have used knockout mutants of each gene and analyzed the mutational effects on the filament length and on the composition and localization of each flagellin species in the filament by electron microscopy and one- and two-dimensional polyacrylamide gel electrophoresis. The location and amounts of flagellins in a filament were determined to be as follows: a small amount of FliC3 at the proximal end, followed by a large amount of FliC5, a large amount of FliC1, a small amount of FliC2 in this order, and a large amount of FliC6 at the distal end. FliC4 was present at a low level, but the location was not determined. Filament lengths of newly born progeny cells increased during prolonged incubation in nutrient-deficient buffer. The newly formed part of the elongated filament was composed of mainly FliC6. Reverse transcription PCR analysis of flagellar gene expression over 5 days in buffer showed that fliC gene expression tailed off over 5 days in the wild-type cells, but in the fliC5 mutant, expression of the fliC2, fliC4, and fliC6 genes was elevated on day 5, suggesting that they may be expressed to compensate for the absence of a major component, FliC5.

Comparative analysis of Klebsiella pneumoniae genomes identifies a phospholipase D family protein as a novel virulence factor

BACKGROUND: Klebsiella pneumoniae strains are pathogenic to animals and humans, in which they are both a frequent cause of nosocomial infections and a re-emerging cause of severe community-acquired infections. K. pneumoniae isolates of the capsular serotype K2 are among the most virulent. In order to identify novel putative virulence factors that may account for the severity of K2 infections, the genome sequence of the K2 reference strain Kp52.145 was determined and compared to two K1 and K2 strains of low virulence and to the reference strains MGH 78578 and NTUH-K2044.

RESULTS: In addition to diverse functions related to host colonization and virulence encoded in genomic regions common to the four strains, four genomic islands specific for Kp52.145 were identified. These regions encoded genes for the synthesis of colibactin toxin, a putative cytotoxin outer membrane protein, secretion systems, nucleases and eukaryotic-like proteins. In addition, an insertion within a type VI secretion system locus included sel1 domain containing proteins and a phospholipase D family protein (PLD1). The pld1 mutant was avirulent in a pneumonia model in mouse. The pld1 mRNA was expressed in vivo and the pld1 gene was associated with K. pneumoniae isolates from severe infections. Analysis of lipid composition of a defective E. coli strain complemented with pld1 suggests an involvement of PLD1 in cardiolipin metabolism.

CONCLUSIONS: Determination of the complete genome of the K2 reference strain identified several genomic islands comprising putative elements of pathogenicity. The role of PLD1 in pathogenesis was demonstrated for the first time and suggests that lipid metabolism is a novel virulence mechanism of K. pneumoniae.

Proteomics and in silico approaches to extend understanding of the glutathione transferase superfamily of the tropical liver fluke Fasciola gigantica

Fasciolosis is an important foodborne, zoonotic disease of livestock and humans, with global annual health and economic losses estimated at several billion US$. Fasciola hepatica is the major species in temperate regions, while F. gigantica dominates in the tropics. In the absence of commercially available vaccines to control fasciolosis, increasing reports of resistance to current chemotherapeutic strategies and the spread of fasciolosis into new areas, new functional genomics approaches are being used to identify potential new drug targets and vaccine candidates. The glutathione transferase (GST) superfamily is both a candidate drug and vaccine target. This study reports the identification of a putatively novel Sigma class GST, present in a water-soluble cytosol extract from the tropical liver fluke F. gigantica. The GST was cloned and expressed as an enzymically active recombinant protein. This GST shares a greater identity with the human schistosomiasis GST vaccine currently at Phase II clinical trials than previously discovered F. gigantica GSTs, stimulating interest in its immuno-protective properties. In addition, in silico analysis of the GST superfamily of both F. gigantica and F. hepatica has revealed an additional Mu class GST, Omega class GSTs, and for the first time, a Zeta class member.

Functional characterization of SjB10, an intracellular serpin from Schistosoma japonicum

Serine protease inhibitors (serpin) play essential roles in many organisms. Mammalian serpins regulate the blood coagulation, fibrinolysis, inflammation and complement activation pathways. In parasitic helminths, serpins are less well characterized, but may also be involved in evasion of the host immune response. In this study, a Schistosoma japonicum serpin (SjB10), containing a 1212 bp open reading frame (ORF), was cloned, expressed and functionally characterized. Sequence analysis, comparative modelling and structural-based alignment revealed that SjB10 contains the essential structural motifs and consensus secondary structures of inhibitory serpins. Transcriptional profiling demonstrated that SjB10 is expressed in adult males, schistosomula and eggs but particularly in the cercariae, suggesting a possible role in cercarial penetration of mammalian host skin. Recombinant SjB10 (rSjB10) inhibited pancreatic elastase (PE) in a dose-dependent manner. rSjB10 was recognized strongly by experimentally infected rat sera indicating that native SjB10 is released into host tissue and induces an immune response. By immunochemistry, SjB10 localized in the S. japonicum adult foregut and extra-embryonic layer of the egg. This study provides a comprehensive demonstration of sequence and structural-based analysis of a functional S. japonicum serpin. Furthermore, our findings suggest that SjB10 may be associated with important functional roles in S. japonicum particularly in host-parasite interactions.

On the alignment of debris disks and their host stars' rotation axis -implications for spin-orbit misalignment in exoplanetary systems

It has been widely thought that measuring the misalignment angle between the orbital plane of a transiting exoplanet and the spin of its host star was a good discriminator between different migration processes for hot-Jupiters. Specifically, well-aligned hot-Jupiter systems (as measured by the Rossiter-McLaughlin effect) were thought to have formed via migration through interaction with a viscous disc, while misaligned systems were thought to have undergone a more violent dynamical history. These conclusions were based on the assumption that the planet-forming disc was well-aligned with the host star. Recent work by Lai et al. has challenged this assumption, and proposes that the star-disc interaction in the pre-main sequence phase can exert a torque on the star and change its rotation axis angle. We have estimated the stellar rotation axis of a sample of stars which host spatially resolved debris disks. Comparison of our derived stellar rotation axis inclination angles with the geometrically measured debris-disk inclinations shows no evidence for a misalignment between the two.